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Q . Abstract 

We describe a method for quantitative characterization of radio-frequency (RF) 

,^ \ inhomogeneity using NOON states. NOON states are special multiple-quantum co- 

^-^' herences which can be easily prepared in star-topology spin-systems. In this method 

■y ' we exploit the high sensitivity of the NOON states for z-rotations. As a result, Torrey 

C^ \ oscillations with NOON states decay much faster than that of single quantum coher- 

^^1 ences. Therefore the present method requires shorter pulse durations, and enables 

the study of inhomogeneity at higher RF powers. To model such a inhomogeneity 

profile, we propose a two-parameter asymmetric Lorentzian function. The experi- 

^ \ ments are carried out in ^H channels of two different high-resolution NMR probes 

'sj" ■ and the results are compared. Finally, we also extend the NOON state method to 

I . characterize the inhomogeneity correlations between two RF channels. We obtain 

tH- \ the correlation profile between ^H and ^^P channels of an NMR probe by construct- 

^f^ ' ing a five-parameter asymmetric-3D Lorentzian model. 
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1 Introduction 



The strength of NMR over other spectroscopy techniques is in the excel- 
lent control over quantum dynamics pLj. Coherent control of nuclear spins 
is achieved by a calibrated set of radio frequency pulses. Due to its low sen- 
sitivity, NMR needs an ensemble of spatially distributed replicas of a spin 
system. NMR signal is the ensemble average of signals arising from each of 
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these replicas. Usually, it is assumed that each member of the ensemble un- 
dergoes the same unitary dynamics. In practice, however, RF amplitude has 
a spatial distribution over the sample volume. This is usually referred to as 
RF inhomogeneity (RFI) [2]. 



Estimation of RFI in a given NMR probe is important not only to char- 
acterize the quality of the probe, but also for experiments requiring precise 
control of spin-dynamics. RFI characterization has helped designing robust 
and high-fidelity quantum gates for quantum information studies [3]. In MRI, 
RFI characterization can help to understand certain image distortions and to 
correct them HllSl. 



The standard method for the measurement of RFI is via Torrey oscillation 
[6]. A single, on-resonant, low-power RF pulse of variable duration is applied 
on a spin-1/2 system with suitably long relaxation time constants. The pulse 
diagram is shown in Fig. [T]^a). Ideally, the magnetization should nutate, with- 
out decay, about the direction of RF. But because of RF inhomogeneity, the 
magnetizations at different parts of the sample nutate with slightly different 
frequencies, ultimately resulting in the decay of the overall magnetization. 
This decaying spiral in vertical plane perpendicular to the direction of RF 
is known as Torrey oscillation. Since NMR signal is always proportional to 
the transverse component of the magnetization, the amplitude of the signal 
collected at different points of the Torrey oscillation also displays oscillatory 
decay. The Fourier transform of this decay can be used to characterize RFI. 
While this method seems to be the simplest way to characterize RFI, the main 
disadvantage is the necessity of long RF pulses to allow complete decay of Tor- 
rey oscillation. Fourier transform of incomplete decays will lead to artifacts. 
In a good NMR probe, Torrey oscillation at the maximum allowed RF power 
may last longer than several milli-seconds and exceeds the duty-cycle limit. 
This makes the measurement of RFI at high powers prohibitive. 



Here we describe a method to speed-up the decay of Torrey oscillations using 
NOON states. NOON states are the highest order multiple-quantum coher- 
ences in a given spin-system. Previously, NOON states have been used in NMR 
for ultra-low magnetic field sensing [7] and for diffusion studies [8]. 



In the following section we describe the theory of Torrey oscillation, prepa- 
ration of NOON states, the measurement of RFI using NOON states, and 
characterization of RFI correlations between two channels of a high resolution 
NMR probe. The experimental details are described in section 3, and finally 
we conclude in section 4. 



2 Theory 



2. 1 Radio frequency inhomogeneity 



For a given nominal RF amplitude uq, let us assume RFI distribution {i/, p(i/)}, 
where p(z/) is the probability, in the sample volume, of the RF amplitude z/. 
The mean RF amplitude is 17 = Y^,^p{i^)i^- Here we have defined the nominal 
RF amplitude z/q as the one with the highest probability. In an asymmetric 
RFI distribution, uq may be different from z7. 

Consider the standard experiment for estimating RFI via Torrey oscillation 
as shown in Fig. [T]^a). During an on-resonant pulse, the distribution of nu- 
tation frequencies ui of the single quantum coherence is same as the RFI 
distribution p(z/). Therefore real part of the signal collected at the end of a 
variable-duration (r) y-pulse is s(r)e~*/"^2, where 
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Fig. 1. (a) Pulse sequence for single-quantum Torrey oscillation, (b) star-topology 
spin-systems, (c) pulse sequence for NOON Torrey oscillation, and (d) pulse se- 
quence for the measurement of correlation between two RF channels. In (c) and 
(d), except the variable duration (r) pulses, all other pulses are 90° pulses. The 
operators corresponding to the dashed pulses cancel each other and can be omitted. 



s{t) = "^pii^) sin(27rz/r). (1) 

V 

Here T2 is the spin-spin relaxation time constant. Similarly, for a g-quantum 
coherence 



s r 



^p(z/) sin(27rgz/r). (2) 



The Fourier transform along t of s(r)e~*/^2 will result in a Lorentzian of line- 
width 7r/T2 centered at zero frequency (on- resonant case) [9] . The amplitude of 
this Lorentzian, is proportional to s(r). We record a series of experiments by 
systematically incrementing r and obtain the Lorentzian peak- heights sir). 
In a typical NMR probe, the real part of the Fourier transform 

S{y,) = Tr \s{r)\ (3) 

displays an asymmetric profile, with a maximum at the nominal RF frequency 
vq [3]. This is because in most cases, the probability of finding lower-than- 
nominal amplitudes is more compared to higher-than-nominal. This profile of 
RF amplitudes can be modeled by an asymmetric Lorentzian 
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The area of p(z/) is normalized to unity. 

Using such a model for RFI it is possible to reconstruct g-quantum Torrey 
oscillations 



v=2vq 

Sm(r) = ^ p(z/)sin(27rgi/T), (5) 

i/=0 

where v is the RF frequency with probability "piy). The positive real part of the 
Fourier transform of Sm(r) gives the model profile S^iyq}- A scalar measure 
of distance between the model and the experimental profiles is provided by 
the norm ||5'(z/q) — Sra{yq)\\- The RFI parameters (Ai, A2) can be obtained by 
minimzing the above norm using a non-linear fit. 

For the hypothetical case of Ai = A2 := A, one can model the Torrey oscil- 
lation as T^^ \Sm{yq)\ oc exp(— Az/ogr). Therefore increase in the nominal RF 
frequency uq will correspondingly speed-up the decay of Torrey oscillation. 



Another way to speedup the decay of Torrey oscillation is to exploit multiple- 
quantum coherences. In the following we explain how NOON states can be 
used for efficient characterization of RFI. 



2.2 NOON state preparation 



We label the Zeeman eigenstates of a single spin-1/2 nucleus as {|0), |1)}. In 
an A^-spin system, the NOON state is the superposition of a state | A^, 0) (with 
A^ spins in |0) and spins in |1)) and the state |0, A^) (with spins in |0) and 
A^ spins in |1)), 



|Ar00Ar) = (|Ar,0) + |0,Ar))/v^ 
= (|00---0) + I11---1 



(6) 



Consider a star-topology spin system as shown in Fig. [T](b). Essentially it is 
an AMi^-i spin system, with a single heteronucleus A interacting with A^ — 1 
magnetically equivalent M spins via J-coupling. An example system is shown 
in Fig. [21 The spectrum of A spin is split into A^ lines, depending on the spin- 
states of M^r-i spins. One of the extreme lines of this spectrum corresponds 
to the coherence 4j(|0) -|- |1))|00 ■ ■ ■ 0), where each M spin is in the state |0) 
(see Fig. 2). This coherence can be converted into a NOON state using CNOT 
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Fig. 2. Molecular structure of trimethylphosphite and its ^^P spectra in equilib- 
rium (upper trace) and after converting NOON state into single quantum coherence 
(lower trace). 



gates {NOT on all M spins controlled by A spin), 



i=(|0) + |l))|00---0)^-^^i=(|00---0) + |ll---l)) = |iVOOiV). (7) 

Since NMR detection is confined to single quantum coherences, it is necessary 
to convert the NOON state back to single quantum coherence before detection. 
This can be achieved by the same set of CNOT gates 



|Ar00Ar)^-^^^(|0) + |l))|00---0). (8) 

v2 

The pulse diagram for NOON state preparation and selection is shown in Fig. 
Wi,c). It can be observed that the CNOT gates in star-topology systems can be 
implemeted in parallel, and therefore require just two RF pulses on M spins. 
A pair of PFGs, Gi and G2, serve to select the above coherence pathway, and 
extract the corresponding extreme spectral line (see Figj2]). The PFG ratio 



^ = (iV-l)^ + l, (9) 

depends on the size of the spin system (A^) and the relative gyromagnetic ratio 

7m /7a- 



2.3 Measurement of RFI via NOON states 



The pulse diagram for the measurement of RFI via NOON states is shown 
in Fig. Wic). In order to measure RFI without changing the coherence order 
of NOON states, we utilize z-pulses. A (j)^ pulse acting on a single spin-1/2 
system introduces a relative phase shift 



;^(|0) + |l))^^(|0) + e^*^|l)). (10) 

If the same pulse is applied on the M^-i spins of a NOON state, and subse- 
quently converted into single quantum coherence, we obtain 



^(|00---0) + |ll---l)) ^ ^(|00---0) + e(^"i)''^^|ll---l)) 



^-^^^(|0) + e'(^-^)'^^|l))|00---0). (11) 

V2 



Therefore the resultant phase-shift of the extreme spectral line is (p^ = {N — 

1)02- A 0^ pulse, on M-spins can be realized by f|j 0^^ f | j . Here 0*^ = 

27t{N — l)i'oT is the corresponding x-pulse of duration r for a nominal RF 
amplitude z/q. Ideally, in the absence of RFI, one should see a regular oscillatory 
behavior of this coherence with r. In practice we see Torrey oscillation, i.e., 
decaying oscillations, due to the RFI during 0;^ pulse. As described in section 
2.1, NOON Torrey-oscillations in large spin-systems decay much faster than 
the single-quantum Torrey oscillations. The faster decay of NOON Torrey- 
oscillations allows one to study RFI at higher amplitudes. 



Here RFI of other short and fixed-duration pulses only introduce an over-all 
scaling factor, and do not contribute to the decay constant of signal. Also, we 
assumed that the transverse relaxation time constant of NOON state is much 
longer than the decay constant of NOON Torrey oscillations. 



2.4 Correlation between RFI of two channels 



In a two-channel probe, the regions of high RF intensity of first channel may 
not correspond to regions of high RF intensity of the second. In other words, 
there exists certain correlations between the RFI profiles of the two chan- 
nels. The NOON state method can be easily extended to study such corre- 
lations in the inhomogeneities of two different RF channels. We choose ^H 
and ^^P channels of Bruker QXI probe for studying such correlations and use 
trimethylphosphite (see Fig. [2]) as the star-topology spin system {^^P^Hg). 
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Fig. 3. Convention used for selection of /3 parameters of the 3D asymmetric 
Lorentzian. The central point (uq,!^^) indicates the nominal point, and the four 
quadrants have different pairs of /3 parameters. 



The pulse-diagram for the RFI-correlation study is shown in Fig. [Hd). A 0f 
pulse is introduced after the (j)^ pulse, and the two pulses are independently 
incremented to obtain a 3D dataset. Similar to Eq. [31 the Fourier transform 
along the two dimensions result in signal S{i>g , i/[). We now model the RFI 
correlation using a 3D asymmetric Lorentzian: 



P(^^,^n= ...... .°p. , ,. , (12) 



where d{u^ , u^) is the scaled distance of ( ^, ^ ) from the nominal point 

(1,1): 
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Here f3^ , /9|^, /3f, and /3^ are the four asymmetry parameters, on the four 
quadrants of z/^ - u^ plane (Fig. [3]). The probability distribution is normalized 
to unit sum. 

As described in section 2.1, we first constrcut the model torrey oscillation 
SmiT^ , T^) and its Fourier transform Sm{i^g , i^i)- By fitting Smi^^g , ^i) with 
the experimental profile S^Uq , z/f), we obtain the best fit parameters (Aq, {/3})- 

In the following, we describe experimental characterizations of (i) ^H RFI at 
different RF powers in two NMR probes, and (ii) RFI correlation between ^H 
and ^^P channels of a NMR probe. 



3 Experiments 



The single- quantum Torrey oscillations were studied using a sample consisting 
of 600 ^1 of 99% D2O. The NOON Torrey oscillations were studied using 100 
|j.l of trimethylphosphite (Fig. [2]) dissolved in 500 \x\ of dimethylsulphoxide- 
D6. The J-couplings between the magnetically equivalent methyl protons and 
the ^^P spin are all equal to 11 Hz. The effective Tg relaxation time con- 
stants of ^H and ^^P are about 150 ms and 730 ms respectively (as is the 
general case. Tip is longer than Tg). These values were obtained from respec- 
tive spectral linewidths. The NOON state T^ is about 38 ms, and is obtained 
by introducing a variable delay between the two CNOTs, and then measuring 
the decay constant of the single quantum coherence after conversion. All the 
experiments were carried out on a Bruker 500 MHz NMR spectrometer at an 
ambient temperature of 300 K. 
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Fig. 4. Characterization of RFI in ^H channel of BBO probe. The results with 
single quantum method are in the top trace, and other traces are for NOON states 
with different nominal RF amplitudes (as indicated in column (a)). Column (a): 
Torrey oscillations s(t) vs t. Column (b): Positive real part of Fourier transform 
S{i'^) vs u^ (solid lines). Also shown are the best fits with model Sm{i^q) (dashed 
lines). Column (c): The RFI profiles p(i^ /i^g ) vs u /vq . The asymmetric linewidth 
parameters (Ai,A2) are indicated in each case. 
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Fig. 5. Characterization of RFI in ^H channel of QXI probe. The results with 
single quantum method are in the top trace, and other traces are for NOON states 
with different nominal RF amplitudes (as indicated in column (a)). Column (a): 
Torrey oscillations s(t) vs t. Column (b): Positive real part of Fourier transform 
S{v^) vs u^ (solid lines). Also shown are the best fits with model Sm{vq) (dashed 
lines). Column (c): The RFI profiles ^{i'^ /uq) vs u^ /uq . The asymmetric linewidth 
parameters (Ai,A2) are indicated in each case. 
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3.1 RFI characterization of ^H channel 



The single quantum Torrey oscillations were studied using the pulse sequence 
shown in Fig. 1(a). A series of spectral lines were recorded by incrementing the 
duration of the low-power (about 1 kHz) on-resonant pulse on 99% D2O sample 
as described in section 2.1. The results of the single quantum experiments on 
^H channels of BBO and QXI probes are shown in the top trace of Figs, 
m and \5\ respectively. The nominal RF frequencies were found to be slightly 
different for the two probes as indicated in the figures. The r-increments in 
BBO and QXI were set to 112.5 [is and 250 |js respectively, and a total of 256 
transients were recorded. The Torrey oscillation (s(r)), the positive real part 
of the Fourier transform (S'(i/^)) are shown in the first traces of Figs. Ill^a,b) 
and[SI^a,b). By fitting the model Sm{i^i) to S^u^), we obtain the asymmetry 
parameters Ai and A2 and the corresponding RFI profile p{v^). The top traces 
of Figs, m^c) and [5]^c) show the RFI profiles p{v^) and the corresponding A 
values for BBO and QXI probes respectively. 



The NOON state is prepared with the ^^P-^Hg system of trimethylphosphite 
(Fig. [2]) using the pulse sequence shown in Fig. Wis). The results of the NOON 
Torrey oscillations at different nominal amplitudes of ^H channels of BBO 
and QXI probes are shown in Figs. Hl^a) and[5]^a). It can be observed that as 
the nominal amplitude is increased, the Torrey oscillations decay faster. The 
positive real parts of the Fourier transform S{v^) of the 9-quantum Torrey 
oscillations are shown in Figs. 111(b) and[5]^b). By fitting the model function 
Smi^g) to S^Ug) at each of these nominal amplitudes, we obtained the cor- 
responding RFI profiles p{v^). These profiles and the asymmetry parameters 
(Ai, A2) are shown in Figs. |l](c) and[S](c). 



It can be observed that decay times of the NOON state Torrey oscillations are 
an order of magnitude shorter than the single quantum Torrey oscillations in 
both probes. This fact allowed us to characterize RFI safely even at maximum 
allowed RF powers. The excellent fit in all the cases indicates the validity of 
our model. In the overall comparison, we find that the BBO probe has stronger 
RFI than QXI probe at all RF amplitudes. The mean values of (Ai,A2) are 
(0.057,0.012) and (0.019,0.009) respectively in BBO and QXI probes. 



In QXI probe, the RFI profiles are not only sharper than those in BBO, but 
are also more symmetric. The mean ratios A1/A2 are 5.2 and 2.2 in BBO and 
QXI respectively. QXI probe shows some what stronger RFI at higher RF 
amplitudes. For example, RFI profiles at 9.5 kHz and 21 kHz are significantly 
wider than the profile at 5.4 kHz. 
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3.2 RFI correlation between ^H and ^^P channels 



We measured RFI correlation between ^H and '^^P channels of Bruker QXI 
probe using the NOON state method. The pulse sequence for measuring 
RFI correlation between ^H and ^^P channels is shown in Fig. [U^d) and 
the theory is described in section 2.4. The dataset of 2D Torrey oscilla- 
tions was recorded by independently incrementing the t^ and r^ pulses re- 
spectively by 11.1 \xs and 50.0 \xs. The nominal RF amplitudes in ^H and 
^^P channels were 2.4 kHz and 2.5 kHz respectively. A total of 128 data 
points in ^H dimension and 96 data points in ^^P dimension were collected. 
The real positive part of the Fourier transform of the 2D Torrey oscillations 
S{v^ , z/f ) was obtained after a zero-fill to 256 points in each dimension. The 
contour plot of S{v}^ .,v^) is shown in Fig. [!](&)■ The asymmetric Lorentzian 
parameters (Aq, {/3}) were obtained by fitting the model profile Sm{u^ .u^) 
to ^(//(f , z/f ). Fig. [6]^b) displays Smiy^ juf) corresponding to the best fit: 
{Ao = 0.0045, /3f = 0.114, /3f = 0.226, /3f = 0.095, /3f = 0.028}. It can be 
seen that the model profile Sm{i-'^ ii^i) approximates the overall profile of 
S^u^^v^). The corresponding RFI profile p{v^ ,v^) is shown in Fig. El^c). 
It can be observed that ^^P-channel has a wider RFI distribution than the 
^H-channel. 



4 Conclusions 



We have proposed an asymmetric Lorentzian model to describe radio- frequency 
inhomogeneity. Using this model, we obtained the asymmetric line-width pa- 
rameters in two high-resolution NMR probes using single quantum Torrey 
oscillations. The excellent fit between the experimental profile and the simu- 
lated profile indicated the suitability of the asymmetric Lorentzian model. 

We have described the characterization of RFI using NOON states. NOON 
state Torrey oscillations decay much faster than single quantum Torrey oscil- 
lations allowing the characterization of RFI at higher RF amplitudes. Using 
this method, we have studied RFI of two NMR probes at different RF ampli- 
tudes and compared the results. We have also proposed a 3D Lorentzian model 
to describe RFI correlations between two RF channls. Then we extended the 
NOON state method to characterize such correlations, and demonstrated it 
on ^H-'^^P channels. We believe that these characterizations will be useful to 
design high-precision NMR and MRI experiments as well as to design high- 
fidelity quantum gates required in NMR quantum information processing. 
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Fig. 6. Characterization of RFI correlations between ^H and '^^P channels of QXI 
probe, (a) Contour plot showing real positive part of Fourier transform of 2D Torrey 
oscillations S{v^ ,vi), (b) corresponding simulated profile obtained with the best 
fit model Sm{i^Q , z^f ), (c) the surface plot of the RFI profile p{v^ ■, v^), and (d) the 
contour plot oi p{v ,v ). 
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